Interaction of a two-level atom with a classical field in the context 

of Bohmian mechanics 



S. V. Mousavi* and M. Golshanit 
Department of Physics, Sharif University of Technology, 
P. 0. Box 11365-9161, Tehran, Iran and 
Institute for Studies in Theoretical Physics and Mathematics (IPM), 
P. 0. Box 19395-5531, Tehran, Iran 

Abstract 

We discuss Bohmian paths of the two-level atoms moving in a waveguide through an external 
resonance-producing field, perpendicular to the waveguide, and localized in a region of finite diam- 
eter. The time spent by a particle in a potential region is not well-defined in the standard quantum 
mechanics, but it is well-defined in the Bohmian mechanics. Bohm's theory is used for calculating 
the average time spent by a transmitted particle inside the field region and the arrival-time distri- 
butions at the edges of the field region. Using the Runge-Kutta method for the integration of the 
guidance law, some Bohmian trajectories were also calculated. Numerical results are presented for 
the special case of a Gaussian wave packet. 
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I. INTRODUCTION 



We study the effect of a classical field on the motion of the ultracold two-level atoms in 
the context of Bohmian mechanics, taking into account reflection effects, and based on the 
quantised longitudinal motion. In general, we have to use at least two counter propagating 
fields to avoid the spatial separation of the internal states due to transverse momentum 
transfer. However, we follow [1] to direct the atoms in a narrow waveguide and work in 
a regime for which the excitation of the transversal modes may be neglected. The field 
changes the longitudinal atomic momentum as it acts as a quantum mechanical potential. 
The incident wave packet is divided into two parts once it reaches the field: the reflected 
one and the transmitted one. We consider the passage of the particles through the field. In 
the context of Bohmian mechanics [2-9], the velocity of a particle is given by v = -, where 
J = V ^ — (V^^)^j i s the current density and p = is the probability density. 

The products on the right hand side of the velocity formula are spinor inner products for 
the two-component wave function. 

For a one-dimensional scattering problem in which a particle is incident normally (from 
the left) on the potential U(r) = V(x)Q(x)Q(£ — x), three characteristic times are defined: 
the transmission time tt(cl, b) and the reflection time tr(cl, b) are the average times spent 
in the region a < x < b by transmitted and reflected particles respectively; the dwelling 
time Tn(a, b) is the average time spent between a and b irrespective of whether the incident 
particle is ultimately transmitted or reflected [10]. 

For a brief discussion on the importance of time and time keeping see [11]. There is an 
ongoing debate on defining the time that characterizes the passage of a quantum particle 
through a given region. For a recent review see [12]. In quantum mechanics, as opposed 
to classical physics, the meaning of the arrival-time of a particle at a given location is not 
evident, when the finite extent of the wave function and its spreading becomes relevant, 
as is the case for cooled atoms dropping out of a trap. Moreover, in the quantum case 
one expects an arrival-time distribution, and there are different theoretical proposals for 
it (see e.g. the review article in [13], and the book [14].) As pointed out by Hannstein 
et al [15] these arrival-time distributions have been controversial, since they are derived by 
purely theoretical arguments without specifying a measurement procedure. They provide an 
operational distribution of arrival times obtained by fluorescence measurements and compare 
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or relate it with more ideal or abstract results. 

Within Bohmian mechanics, the average characteristic times are uniquely defined and 
each one is obviously a real, non- negative and additive quantity [16]. As Bohm has empha- 
sized in connection with the measurement of the momentum of a particle [2-4] , the quantity 
actually measured may have no relation to its value before the measurement. In Leavens's 
view [16] the basic difference between the Bohm's trajectory approach and the 'conventional' 
approach to characteristic times is that Bohmian mechanics consistently treats a 'particle' 
as an actual particle at all times, not just at those instants when an ideal (i.e. strong) 
measurement is made of its position. 

As the numerical calculation involved in the calculation of Bohmian trajectories is im- 
mense, there have been attempts to compute the Bohmian transmission and reflection times 
without sampling the trajectories [17-19]. 

The next section contains the solution of the stationary Schrodinger equation. This 
is followed by a brief review of the relevant parts of Bohm's interpretation of quantum 
mechanics. Sec. IV, contains some numerical calculations. Concluding remarks are made in 
Sec. V. 

II. SOLUTION OF THE STATIONARY SCHRODINGER EQUATION 

The most common interaction between an atom and an electromagnetic field is of the 
electric dipole form which couples the states |1) and |2) of differing parity. The electric dipole 
approximation is valid when the wavelengths of the driving field are much greater than the 
mean separation of the electron and the nucleus. As quantum optics is primarily concerned 
with the interaction between atoms and optical or infrared radiation, this approximation is 
usually a good one, because the approximation conditions can be met in the laboratory. The 
atomic dipole operator associated with the two states |1) and |2) has the general form of /i = 
/x*|2)(l| + /i|l)(2|, where the caret^is used to distinguish operators from the corresponding 
c-numbers. The Hamiltonian describing a two-state atom interacting with a classical electric 
field E(t), within the electric dipole approximation, is [20] 

H = £ + |k|1)(1|+c 2 |2)(2|)-/2.E, (1) 
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where we have taken /i to be real for simplicity. The quantities Huj\ and huji are the energies 
of the states |1) and |2) respectively. In the interaction picture, and using the rotating-wave 
approximation, Hamitonian takes the form [1] 



H = ^--hA\2)(2\ + ^n(x)(\2)(l\e-^ + h.c), (2) 

where Q(x) is the position-dependent Rabi frequency, is the phase of the classical field, and 
A = ujl — oj\2 denotes the detuning between the field frequency and the atomic transition 
frequency, u\2 = —Wi, and h.c denotes hermite conjugate. For simplicity we assume that 
Q(x) = Q for < x < £, and zero otherwise. To obtain the time development of a wave 
packet incident from the left (k > 0), we first solve the stationary equation HQ^ = E^^ for 
the scattering states with the energy E = 4^ = E k , 



1 / e ikx + R 1 (k)e- ikx , 
= = -?= I \ x < 



2vr \ R 2 (k)e- iqx 



1 / Ti(k)e ikx 

= ®kA x ) = I x > £ ■ ( 3 ) 

V27T I T 2 {k)e iqx 



The wavenumber q obeys q — y k 2 + due to the conservation of energy, while R\^{k) 
and Ti^ (fc) are reflection and transmission amplitudes for the energy E k . Within the field 
region, the (unnormalized) dressed state basis that diagonalise the Hamiltonian is given 
by |A±) = |1) + ^j^e~ l< ^\2) , where X± = (— A±fi')/2 are the dressed eigenvalues in which 



VI' = y/ A 2 + Vt 2 . Thus, the solutions in the interaction region is of the form, 



|$ fc (a0) = -j= ^A {+) e ik+x + B {+) e- ik+x }\\ + ) + (A^e ik ~ x + B^e - '*-*) |A_) , (4) 

with the wavenumbers k 2 ± = k 2 — 2m\±/K. The continuity of $/c(x) and its derivative at 
x = and x = i leads to eight equations with eight unknowns. To impose the matching 
conditions, it is convenient to use a two-channel transfer matrix formalism [21]. 

For stationary states, the density probability p = $|,(x)<I>fc(x) is independent of time, and 
from the continuity equation in 1-dimension, §f + §i£ = 0, it follows that the current density 
is a constant. For the left and the right of the field region, the current density is: 



■h = -{l-\R 1 {k)?)--\R2{k)\\ 
m m 

J R = ^{k)? + ~\T2{k)\\ (5) 
m m 



and from Jl = Jr one has 

|-Ri(^)| 2 + l^k)] 2 + ^(\R 2 (k)\ 2 + \T 2 {k)\ 2 ) = 1 . (6) 
The wave packet which describes the particle is of the form 

q(x,t) = J dk^ k (x)^(k)e- iEkt/h , (7) 

where we choose ip(k) to be 

#( k ) = (^0)V4 e -^(fc-fe ) 2 e -to Q ; (g) 

7T 

in which Hka is the momentum kick as the result of which the particle moves along x— axis 
towards the field, x is the centre of the wave packet at t — 0, and a is the initial width of 
the wave packet. 



III. SOJOURN TIMES IN BOHMIAN MECHANICS 

In nonrelativistic Bohmian mechanics the world is described by point-like particles which 
follow trajectories determined by a law of motion. The evolution of the positions of these 
particles is guided by a wave function which itself evolves according to the Schrodinger's 
equation. Bohmian mechanics makes the same predictions as the ordinary nonrelativis- 
tic quantum mechanics for results of any experiment, provided that we assume a ran- 
dom distribution for the initial configuration of the system and the apparatus in the 
from of p(x,0) = ^(x, 0)*f?(x, 0). If the probability density for the configuration satis- 
fies p(x,to) = ^(x, t )ty(x, t ) at some time to, then the density to which this is carried by 
the continuity equation at any time t is also given by p(x,t) = &(x,t)ty(x,t) [22]. Because 
most quantum measurements boil down to position measurements, Bohm's theory and the 
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standard quantum mechanics will in general yield the same detection probabilities. The 
situation is different, however, if one considers, for example, measurements involving time 
related quantities, such as time of arrival, tunnelling times etc. Bohm's theory makes un- 
ambiguous predictions for such measurements, but in the conventional quantum mechanics 
there is no consensus about what these quantities should be [23]. Given the initial position 
= x(t = 0) of a particle with the initial wave function ^(x, t = 0), its subsequent trajec- 
tory x(x^°\ t) is uniquely determined by the simultaneous integration of the time dependent 
Schrodinger equation, TDSE, and the guidance equation = v(x(t),t). 

As pointed out by Leavens, the noncrossing property of Bohmian trajectories x(x^°\t) 
leads to a very important consequence in 1-dimension: there is a critical trajectory (starting 
at Xc) that divides the wave packet into two parts, the transmitted one (x^ > Xc ) and 
the reflected one (x^ < Xc) [24]. The reflection-transmission bifurcation curve x c (t) 
separating the transmitted trajectories from the reflected ones is given by 

POO 

\T\ 2 = / dx¥(x,t)V(x,t) , (9) 

Jx c (t) 

/oo 
dx^' (x, t)^(x, t) and 

could be calculated from the stationary-state transmission probabilities \Ti(k)\ 2 and |T 2 (fc)| 2 
as 



2a, 



2 poo 




dk 

) 

x (|T!(A;)| 2 + ||T 2 (A;)| 2 )e- 2,T o( fc - fc o) 2 , (10) 
The dwelling time td(0,£) in the field region is given by 



poo /•£ 

dt dx^(x,t)^(x,t) . (11) 
Jo Jo 



Inserting 1 = 6(x — x c (t)) + 0(x c (t) — x) into the integrand of this equation, one gets 



r D (0,£) = \T\ 2 tt{Q,£) + \R\ 2 t r (0J) , (12) 
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in which 



r T (0,£) 



T 



dt / dx&(x,t)V(x,t)Q( 
Jo 



x - x c (t)) 



(13) 



poo p£ 

r R (0,£) = jj^J & J dxV\x,t)V(x,t)Q(x e (t) - x) 



(14) 



where td(0, £), tt(0,£) and r R (0,£) are, respectively, the mean dwelling, transmission and 
reflection times respectively, and |T| 2 and |i?| 2 = 1 — |T| 2 are transmission and reflection 
probabilities, respectively. Within the Bohmian mechanics, the average characteristic times 
tx=d,t,r(0, £) are uniquely defined and each one is obviously a real, non-negative and additive 
quantity in the sense of measures: Tx(xi,x 3 ) = Tx(xi,x 2 ) + Tx(x 2 ,x 3 ), with xi < x 2 < x 3 , 
in which X could be D, T or R [16]. 

The total presence probability at the right of the point x at time t, Q(x,t), is given by 



poo rt 

Q(x,t) = / dx'¥(x',t)^(x',t) = / dt'J{x,t') 

Jx JO 



(15) 



As discussed by Oriols, Martin and Sune [18] , the problem of the average characteristic times 
can be reduced to the calculation of the total presence probability at the right of the edges 
of the field region, Q(0, t) and Q(£, t): 



I dt[Q(0,t) - Q(£,t) 



(16) 



and 



\T\ 



dt 



min[Q(Q,t), \T\ 2 ) - Q(£,t) 



(17) 



Tr 



\R\ 



dt 



max[Q(0,t), \T\ 2 ) - \T 



(18) 
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In the causal interpretation, the distribution of arrival-times Pb(t) at a given interface 
x = b, with b > £, is related to the current density J(b, t) by 



P(f ) - J{h ^ - 1 Wilt) 

provided there are no re-entrant trajectories through x = b [17, 24]. This expression is 
valid in the absence of an arrival-time detector, because an arrival-time detector is not 
included in the Hamiltonian to derive this expression. Thus, this is an expression for the 
ideal or intrinsic arrival-time distribution [25]. If Q(b,t) is not monotonic, then J(b,t) must 
be negative for some times, implying re-entrant trajectories through b. In this case the 
arrival-time distribution is given by 



m) = f jJ(M)l _ (20) 
dt\J(b,t)\ 



allowing for multiple arrival-times, both from the left and right, associated with any re- 
entrant trajectory that might occur [16-18, 25]. 



IV. NUMERICAL RESULTS 

Consider Cs atoms with mass m = 2.2 x 10 _25 kg, impinging on a classical on-resonance 
|— pulse, with the Rabi-frequency f2 = (n+l/2)x^, where n is an integer number. The 
initial velocity of the particle is chosen to be vq = O.Olm/s, the width and the phase of the 
field is chosen as I = 10 _4 m and = respectively, and n = 650 and oo = 20yum. With 
these values, the fraction of the interaction energy to the kinetic energy of the particle is 
^ = 1.96. 

Fig. 1 shows the probability density for a wave packet approaching the field region and col- 
liding with the field, and also the construction of the transmitted and reflected wave packets. 
This figure resembles the behaviour of a potential well or a potential barrier (with a height 
smaller than the kinetic energy of the arriving particle), when a wave packet approaches the 
potential region. 

Numerical calculations show that transmission probability |T| 2 increases with both do and 
A (Fig. 2 and Fig. 3). Increasing A means that detuning is increased, thus the interaction 
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between the particle and the field reduces, i.e., transmission probability is increased. In the 
case that total energy of the atom, the kinetic energy of the excited state and the interaction 
energy Ml are negligible with respect to HA, one can show, using the stationary Schrodinger 
equation, that the atom propagates freely in the ground state. 

In Fig. 4, a selection of Bohmian trajectories in the neighbourhood of the transmission- 
reflection bifurcation is shown, using the Runge-Kutta method for the integration of the 
guidance equation. From Eq. (10), we get |T| 2 = 0.61 for <jq = 20/im. Since the initial wave 
packet is an even function of x around Xq, the centre of the initial wave packet, the total 
probability at the right and left of xq is each one equal to 0.5. Thus, Eq. (9) at t = and 
\T\ 2 = 0.61 yields Xc — Xq — 0.56 x 10 05 < x , i.e., the bifurcation trajectory starts in the 
back half of the initial wave packet. At points where the effective force (the classical and the 
quantum ones) is small, or equivalently probability density is large, trajectories congregate 
and give rise to fringes, a matter- wave analogue of Wigner's optical fringes [5]. It should be 
noted that the continuity of the wave-function and its derivative at the boundaries, imply 
that probability density and the probability current density both are continuous functions at 
the boundaries. Thus, it is found from the guidance equation that the velocity of a Bohmian 
particle is continuous at the boundaries. This point is not evident from Fig. 4 because of its 
scale, but it was found by zooming on this figure at the boundaries. 

Using Eqs. (16, 17 and 18), one gets r D = 40 ms, tt = 39 ms and t r = 42 ms for cr = 20 
micron. These values satisfy the weighting relation td = |T| 2 rx + \R\ 2 tr, which is an 
unprovable relation in the standard quantum mechanic. Eq. (12) shows that the integrated 
density under the field region consists of additive contributions from the transmitted and 
from the reflected beams. As pointed out by Landauer and Martin [26], one does not sum 
probabilities in quantum mechanics; rather, it is the complex amplitudes that are summed 
over. Eq. (12) neglects the possibility of interference between the amplitudes for reflection 
and transmission. 

Fig. 5 shows the current and the logarithm of the current densities and Fig. 6 shows 
the corresponding integrated current densities, Q(0,t) and Q(£,t), at the edges x = and 
x = as a function of time. We have plotted the logarithm of the current densities to 
see the number of peaks of the current densities. Fig. 5 shows that J(0, t) is negative at 
some instants, corresponding to some re-entrant trajectories at this boundary, but J(£, t) is 
positive for all the times. 
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Finally, Fig. 7 shows the arrival-time distributions at x = and x — £, according to 
Eq. (19) and Eq. (20) respectively. 

V. DISCUSSION 

The calculations presented here in the context of Bohmian mechanics, are for the traversal 
time, the dwelling time and the arrival-time distributions in the absence of any measuring 
device. Hence, before comparison with experiment can be seriously considered, a method 
for doing the measurement must be specified and the effect of the interaction between the 
propagating particle and the apparatus must be included in the calculation. As Bohm 
emphasized, the quantity actually measured may have no relation to the quantity before 
the measurement. An important feature of the causal interpretation is that, since it gives a 
well-defined, unambiguous and unique prescription for calculating the characteristic times 
and the arrival-times, it is refutable experimentally, provided that the Bohmian trajectory 
calculation of the desired quantity and the corresponding experiment can both be carried out 
with sufficient accuracy. On the other hand, since there is no prescription in the basic tenets 
of the conventional interpretation for calculating the characteristic times, that interpretation 
isn't refutable via comparison of the calculated and the actual experimental characteristic 
times [10]. 
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Figure captions 

Figure 1: Probability density versus distance x(m) at times a) t — ms, b) t — 8.7 ms, 
c) t = 14.7 ms, d) t = 23.7 ms, e) t=29.7 ms and f) t = 104.7 ms. Vertical dashed lines 
show edges of the field. 

Figure 2: a) Transmission and b) reflection probability versus o"o(/xm). 

Figure 3: a) Transmission and b) reflection probability versus A (kHz). 

Figure 4: A selection of Bohmian paths. The position of the field is indicated by the hor- 
izontal dashed lines. The black curve, which starts at x^ = —0.126 mm, is the bifurcation 
trajectory. 

Figure 5: a) J(0,t), b) J(£,t), c) log(| J(0, t)\) and d) log(J(£, t)) versus time(s). From a) 
and c) one can find that in a) there is a big positive peak and then five small negative ones; 
the first represents the wave packet as it enters the field region, the later five represent the 
reflected packets. From b) and d) it's concluded that there are five positive peaks in b). 

Figure 6: Qo(t) (solid line) and Qe(t) (dashed line) versus time(s). The horizontal dotted 
line shows the transmission probability |T| 2 . 

Figure 7: Po(t) (solid line) and Pe(t) (dashed line) versus time(s). 
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FIG. 1: Probability density versus distance cc(m) at times a) t = ms, b) t = 8.7 ms, c) t = 14.7 
ms, d) t = 23.7 ms, e) t=29.7 ms and f ) t = 104.7 ms. Vertical dashed lines show edges of the field. 
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FIG. 2: a) Transmission and b) reflection probability versus <7o(/Uto). 
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FIG. 3: a) Transmission and b) reflection probability versus A (kHz). 
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FIG. 4: A selection of Bohmian paths. The position of the field is indicated by the horizontal 
dashed lines. The black curve, which starts at Xc = —0.126 mm, is the bifurcation trajectory. 
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FIG. 5: a) J(0,t), b) J(£,t), c) log(| J(0, t)\) and d) log(J(£,t)) versus time(s). From a) and c) one 
can find that in a) there is a big positive peak and then five small negative ones; the first represents 
the wave packet as it enters the field region, the later five represent the reflected packets. From b) 
and d) it's concluded that there are five positive peaks in b). 
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FIG. 6: Qo(t) (solid line) and Qe(t) (dashed line) versus time(s). The horizontal dotted line shows 
the transmission probability \T\ 2 . 
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FIG. 7: Po(t) (solid line) and Pe(t) (dashed line) versus time(s) 
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